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Abstract 

The relativistic Faddeev equation for three-nucleon scattering is formulated in momentum space 
and directly solved in terms of momentum vectors without employing a partial wave decom- 
position. Relativistic invariance is achieved by constructing a dynamical unitary representation 
of the Poincare group on the three-nucleon Hilbert space. The exclusive breakup reaction at 
508 MeV is calculated based on a Malfliet-Tjon type of two-body interaction and the cross sec- 
tions are compared to measured cross sections at this energy. We find that the magnitude of the 
relativistic effects can be quite large and depends on the configurations considered. In spite of 
the simple nature of the model interaction, the experimental cross sections are in surprisingly 
good agreement with the predictions of the relativistic calculations. We also find that although 
for specific configurations the multiple scattering series converges rapidly, this is in general not 
the case. 
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Breakup reactions in the proton-dcutcron {pd) system at intermediate energies have 
been studied experimentally quite intensively in recent decades. A prominent set of data 
can be found in the comprehensive overview of the experiments completed at Saturne- 
2 [1]. However, the theoretical interpretation faced and still faces serious challenges. At 
those energies pion production channels arc open and nuclear resonances play a role. In 
contrast to the energy regime below the pion threshold, where high precision nuclcon- 
nuclcon (NN) forces arc established [2l3l4j . and nuclear forces based on effective chiral 
dynamics are being developed [5], the intermediate energy regime, in which the lowest 
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nucleoli resonances play a role, does not yet have forces describing data with compara- 
ble quality [6j. Moreover it can be expected that three-nucleon (3N) forces will play a 
more prominent role in the intermediate energy regime. The tendency for an increasing 
importance of 3N forces with increasing energy has already been seen in 3N observables 
at low energies [7]. 

A consistent treatment of intermediate energy reactions requires a Poincarc symmetric 
quantum theory [S]. In addition, the standard partial wave decomposition, successfully 
applied below the pion-production threshold [9], is no longer an adequate numerical 
scheme due to the proliferation of the number of partial waves. Thus, the intermediate 
energy regime is a new territory for few-body calculations, which waits to be explored. 

The aim of this article is to address two aspects in that list of challenges: exact Poincare 
invariance and calculations using vector variables instead of partial waves. In a series of 
publications [lOj the technique of solving the nonrelativistic momentum-space Faddeev 
equation without partial waves has been mastered, for bound as well as scattering states. 
The Faddeev equation, based on a Poincare invariant mass operator, has been formulated 
in detail in The resulting Faddeev equation has both kinematical and dynamical 
differences with respect to the corresponding nonrelativistic equation. 

The formulation of the theory is given in a representation of Poincare invariant quan- 
tum mechanics where the interactions are invariant with respect to kinematic translations 
and rotations [T^. The model Hilbert space is a three-nucleon Hilbert space (thus not 
allowing for absorptive processes). The method used to introduce the NN interactions in 
the unitary representation of the Poincare group allows to input high-precision NN in- 
teractions }2l3l4j in a way that reproduces the measured two-body observables. However 
in this study we use a simpler interaction consisting of a superposition of an attrac- 
tive and a repulsive Yukawa interaction that supports a bound state with the deuteron 
binding energy. Poincare invariance and S'-matrix cluster properties dictate how the two- 
body interactions must be embedded in the three-body dynamical generators. Scattering 
observables are calculated using Faddeev equations formulated with the mass Casimir 
operator (rest Hamiltonian) constructed from these generators. 

In [llj the driving term in the relativistic Faddeev equation (first order in the two-body 
transition operator) has been used to evaluate pd elastic scattering as well as break-up 
cross sections. This has now been completed by fully solving the relativistic Faddeev 
equation based on the numerical techniques previously used to solve the nonrelativis- 
tic Faddeev equation [TD]. Our calculations converge well up to 2 GeV, indicating the 
applicability of the formulation of the Faddeev equation based on vector variables for 
intermediate energies. We want to point out that the relativistic Faddeev equations with 
realistic spin-dependent interactions have been solved below the pion-production thresh- 
old in |13j a partial wave basis. 

In order to estimate the size of relativistic effects the interactions employed in the 
nonrelativistic and relativistic calculations presented here are chosen to be phase shift 
equivalent. This is achieved in this article by adding the interaction to the square of the 
mass operator |15|16j . In this comparison differences in the relativistic and nonrelativistic 
calculations first appear in the three-body calculations. Those differences are in the choice 
of kinematic variables (Jacobi momenta are constructed using Lorentz boosts rather 
than Galilean boosts) and in the embedding of the two-body interactions in the three- 
body problem, which is a consequence of the non-linear relation between the two and 
three-body mass operators. These differences modify the permutation operators and the 



2 



ofT-shell properties of the kernel of the Faddeev equations. 
The relativistic Faddeev equation is solved in the form 

T(z) = t{z)P + t{z)P{z - Mo)-'T{z), (1) 

where P is a permutation operator and t{z) is the two-body transition operator embedded 
in the three-nucleon Hilbert space defined as the solution to 

t{z)^V + V{z~Moyh{z). (2) 

Here the interaction is given as V" = M — Mq, where M is the three-particle mass oper- 
ator with one two-body interaction alone and Mq the corresponding free mass operator. 
The physical transition amplitudes for elastic scattering and the breakup processes can 
be expressed in terms of T{z). Our aim in this article is to calculate exclusive three- 
body breakup processes in the intermediate energy region and compare relativistic and 
nonrelativistic results. Since we generate the solution of the Faddeev equation from the 
multiple scattering series (resulting from iterating Eq. IIJ), we can also obtain insight into 
the contributions of the various orders of that series. The exact solution to the Faddeev 
equation is constructed from the multiple scattering series using Fade summation. 

Though our two-body force is simple, we want to compare to a ^H(p,2p)n experiment at 
a projectile kinetic energy 508 MeV [lT| to see if our calculation captures essential features 
of the measurement. Differences in the predictions of our relativistic and nonrelativistic 
calculations are already very pronounced at this energy. 

The five-fold laboratory differential cross section for exclusive scattering is given by 

m 

d'ajf ^ ,^.A E{q^)Ed{qo)E{q) 

dnidn2dEi ^^""^ 2klabmd P2i^ - E{pi)) ^ E{p2)iP - pi) ■ P2 

xE{k)^4E^k)+(i^ |(k,q||[/o||¥'d,qo)|'. (3) 

Here t/o = (1 + P)T is the operator for breakup scattering, E is the total energy of the 
system and P its total momentum. 

The results of our calculations are displayed in Figs.[T][31 In Fig. [T] we chose selected 
outgoing proton laboratory angle pairs 6i — 62 from Ref . J7J , which are symmetric around 
the beam axis. The cross sections are plotted against the laboratory kinetic energy of 
one of the outgoing protons. It is interesting to see that for the smaller angle pairs, 
6*1,2 = 38.1° and 41.5°, the relativistic cross sections (solid lines) are considerably larger 
than the nonrelativistic ones (dashed lines). For angle pairs around ^1,2 = 44° and larger 
this reverses and the relativistic cross sections fall below the nonrelativistic ones. A 
similar phenomenon has been observed in Ref. [13] for a projectile energy of 200 MeV, 
where the relativistic Faddeev equation has been solved based on the CD-Bonn NN force. 
Surprisingly, our present relativistic calculations come close to the data [T7]. The fact 
that the first order results nearly coincide with the full calculation may be explained by 
the quasi free scattering (QFS) condition, which however, is realized only for the angle 
pair 9i,2 = 41.5° - 41.4°. 

Configurations in which the outgoing protons are measured asymmetric with respect 
to the beam axis are shown in Fig. [2] Again, the relativistic effects are very pronounced in 
all configurations shown. The angle combination 30.1° — 53.7° is the only QFS condition. 
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In order to illustrate how the cross section is built up by the lowest order terms of the 
multiple scattering series, we show in Fig. [3] the relativistic results for the first order 
in the two-body t-matrix, then add successively one and two rescattering terms, and 
compare to the full calculation. For the QFS condition, 41.5° — 41.4°, rescattering terms 
do not play a role, whereas in the two asymmetric configurations even the 3rd order in 
the multiple scattering series is not quite sufficient to arrive at the full result over the 
entire measured energy range. 

To explore the importance of higher order terms in the multiple scattering series, we 
display in Fig. |4] three configurations for the coplanar star and one for the space star. (In 
the star configurations the interparticle pair angles are 120° in the center of momentum 
(cm.) system). The angles Xq = q ■ qo and Xp = p ■ qo refer to the cm. system with p 
and q being Poincare-Jacobi momenta and qo being the initial particle momentum. The 
angle 4'pq = 0° indicates coplanar star configurations, while (j)pq — 90° indicates a space 
star configuration. For Xq = 1 (one particle being ejected forward along the direction of 
qo) the first order result is completely misleading and about two orders of magnitude 
higher than the fully converged result, which in turn requires more than four orders in 
the multiple scattering series. However, the first order clearly shows the shift in the peak 
position characteristic for relativistic calculations The configuration Xq ~ 0.5 also 
exhibits this shift in the peak position when comparing relativistic and nonrelativistic 
calculation, and again the full result is built up by several rescattering contributions. The 
coplanar star at Xq ~ —1 exhibits the QFS condition (in the laboratory frame particle 1 
would stay at rest), and thus the first order is sufficient. In addition there is no shift in 
the peak position due to relativistic effects, which can be understood as consequence of 
the particle being at rest. In case of the space star, the first order does not contribute, 
and the cross section of this configuration is slowly built up by rescattering contributions, 
requiring orders higher than four. The star configuration shown is representative for all 
space star configurations achieved by rotation in the plane perpendicular to the beam. 

These observations suggest that even at intermediate energies the full solution of the 
Faddeev equation is needed to make reliable predictions. There is no loss of generality 
in using this formalism: because the two-nucleon interactions are fit to two-nucleon data 
and S'-matrix cluster properties are used to embed the two-nucleon interactions in the 
three-nucleon system, any other Poincare invariant three-nucleon model with these prop- 
erties differs from this one by at most a three-body interaction. The difference between 
the relativistic and nonrelativistic calculations using the same two-body input first ap- 
pears in the three-body system. Our calculation indicates measurable differences in the 
corresponding relativistic and nonrelativistic breakup cross sections. Even though our 
model interaction is much simpler than a realistic two-body interaction, the exclusive 
breakup cross section are in surprisingly good agreement with existing data at 508 MeV. 
The results reported in this letter suggest that 

(a) A direct integration without employing a partial wave decomposition can successfully 
be used to solve the relativistic Faddeev equations with embedded NN interactions. 

(b) The formulation of Poincare invariant quantum theory with a three-dimen-sional 
kinematic Euclidean symmetry is a suitable formalism for treating scattering in the 
intermediate energy region. This is supported by the surprising agreement with the mea- 
surements reported in |17| . 

(c) Convergence of the multiple scattering series is not guaranteed at intermediate ener- 
gies. 
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A detailed description of our results comprising elastic, inclusive and exclusive processes 
will follow in a subsequent publication. 
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Fig. 1. The exclusive differential cross section for the ^H(p,2p)n reaction at 508 MeV laboratory projectile 
energy for different proton angle pairs 61-62 symmetric around the beam axis as function of the laboratory 
kinetic energy of one of the outgoing protons. The solid line (R F) represents the full relativistic solution of 
the Faddeev equation, while the dotted curve (R 1) indicates the relativistic calculation based on the 1st 
order in the multiple scattering expansion of the Faddeev amplitude. The corresponding nonrelativistic 
full solution of the Faddeev equation is given by the short-dashed curve (NR F) and its 1st order 
contribution by the double-dotted curve (NR 1). The data are taken from Ref. | 17| . 
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Fig. 2. The exclusive differential cross section for the ■^H(p,2p)n reaction at 508 MeV laboratory pro- 
jectile energy for different proton angle pairs 9i-92 asymmetric around the beam axis as function of the 
laboratory kinetic energy of one of the outgoing protons. The meaning of the curves is the same as in 
Fig. [T] The data are taken from Ref . [T7] . 
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Fig. 3. The exclusive difTcrential cross section for the ^H(p,2p)n reaction at 508 MeV laboratory projectile 
energy for different proton angle pairs 01-02 as function of the laboratory kinetic energy of one of the 
outgoing protons. The solid line (R F) represents the full relativistic solution of the Faddeev equation. 
The dotted line (R 1) represents the relativistic calculation based on the 1st order term in the multiple 
scattering expansion of the Faddeev amplitude, while the short-dashed curve (R 2) includes the first two 
terms, and the dash-dotted curve (R 3) the first three terms. The data are taken from Ref. |17l . 
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Fig. 4. The exclusive differential cross section for the ■^H(p,2p)n reaction at 508 MeV laboratory projectile 
energy for different star configurations as function of the cm. kinetic energy of one of the outgoing protons 
for different cm. angle pairs Xq = cosdq, Xp = cosSp. The angle <j>pq= deg indicates coplanar star 
configurations, while (ppq = 90 deg indicates a space star. All figures show the full relativistic Faddeev 
calculation as solid line (R F). The left column compares the full relativistic calculation with the full 
nonrelativistic one (NR F), as well as the corresponding relativistic (R 1) and nonrelativistic (NR 1) 
first order calculations. The right column shows the subsequent sum of the lowest orders in the multiple 
scattering series contributing to the full relativistic result. The meaning of those curves is the same as 
in Fig. |3] 
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